Identification and characterization of icosahedral metallic nanowires 
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We present and discuss an algorithm to identify and characterize the long icosahedral structures 
(staggered pentagonal nanowires with 1-5-1-5 atomic structure) that appear in Molecular Dynam- 
ics simulations of metallic nanowires of different species subjected to stretching. The use of the 
algorithm allows the identification of pentagonal rings forming the icosahedral structure as well as 
the determination of its number n p , and the maximum length of the pentagonal nanowire L™. The 
algorithm is tested with some ideal structures to show its ability to discriminate between pentago- 
nal rings and other ring structures. We applied the algorithm to Ni nanowires with temperatures 
ranging between 4K and 865K, stretched along the [100] direction. We studied statistically the 
formation of pentagonal nanowires obtaining the distributions of length L™ and number of rings n p 
as function of the temperature. The L™ distribution presents a peaked shape, with peaks locate at 
fixes distances whose separa-tion corresponds to the distance between two consecutive pentagonal 
rings. 

PACS numbers: 02.70.-c, 02.70.Ns, 61.46.Km, 62.23.Hj, 62.25.-g 



INTRODUCTION 

Icosahedral or pentagonal nanowires are formed by 
subsequent staggered parallel pentagonal rings (with a 
relative rotation of 7r/5) connected with single atoms, 
showing a characteristic -5-1-5-1- ordering (see an exam- 
ple in Fig. [TJl Metallic nanowires are of great tech- 
nological importance due to their properties and poten- 
tial applications [UH]- Contrary to monoatomic chains, 
penta-gonal nanowires are rather robust structures at rel- 
atively high temperatures and, therefore, they may con- 
sider as a promising candidate for being used as nanodc- 
vice components. 

Different computational works during the last decade 
have showed the formation of staggered pentagonal con- 
figurations on breaking nanowires of different species 
[3|4l|5j|6j|Il|8j|9j[T0j[ni[l2]- The atomic sequence 
-1-5-1-5- presents a fivefold symmetry with respect the 
nanowire axis. This symmetry does not correspond to 
any crystallographic FCC nor BCC structures. The -1- 
5-1-5- staggered nanowire configuration may be under- 
stood in terms of a sequence of interpenetrated icosahe- 
dra. This icosahedral symmetry is very common in very 
small systems due to the large stability and high coordi- 
nation characterizing such geometry |13j . 

The formation of staggered pentagonal configurations 
during the stretching process has been already reported 
for Na [3] using first principles methods, and for Cu 
[11 [5j |6] and Au [7] nanowires with different Molecular 
Dynamic (MD) approaches. In particular the high stabil- 



ity of the Cu nanowire was confirmed with ab-initio calcu- 
lations [S]. Pentagonal motives also appear in infinite Al 
and Pb nanowires obtained from MD simulated anneal- 
ing methods |8 . More recently such structures have been 
reported for stretched Ni nanowires with different crys- 
tallographic orientations [TUJ [TTJ , and confirmed then- 
stability by ab-initio simulations [HJE2]. These pentag- 
onal structures are very stable, with lengths larger than 
20 A and presenting a high plastic deformation under 
strain. 

In general, the reported pentagonal nanowires have 
been found for single stretching evens. However, it is 
well known that the analysis of nanoscale processes re- 
quires the use of statistical approaches since there exists 
a broad bunch of breaking paths in the nanowire configu- 
ration space. Such statistical studies have been only ad- 
dressed for Ni up to date [TTJl ITTl IT2"] . It has been shown 
that [100] and [110] stretching direction favour the ap- 
pearance of long pentagonal nanowires [lOj [TTJ [12], and 
that there exists an optimal temperature at which the 
pentagonal nanowire yield is maximized [12 . 

In Refs. (TUJ [TTJ [H] a method based on the time that 
the breaking nanowire lasts with a cross section S m ~ 5 
(in units of atoms) was used to detect the formation of -5- 
1-5- structures. The quantity S m ~ 5 is close to the min- 
imum cross-section of a pentagonal ring. As it is shown 
in Fig. [TJd, the formation of the pentagonal nanowire is 
reflected in the curve of the nanowire minimum cross- 
section versus time, presenting a long plateau around 
S m ~ 5 during the pentagonal nanowire formation. A 
large statistical occurrence of pentagonal nanowires (for 



2 



JO Cl SP C 

o> »e © 43 
■ •_« • - 



(a) 



Ni[100] T=375K 



O 63 <(. V I. c '■ ... otf e 

. ° a ° a ° o M-^X S'lTl # f * I ° » 



a 

» 4 4 I 

" • _ « 
_ 9 



"c « c*'c c ci""c u 





600 700 800 900 1000 1100 
t(PS) 



2 4 6 3 10 

S (atomic section) 



FIG. 1: Results from MD simulations of Ni [100] nanowires 
containing 1029 atoms and subjected to longitudinal stretch- 
ing at T=375K. (a) Snapshot of a nanowire formed under 
stretching presenting an icosahedral structure with n p = 10 
pentagonal rings, (b) Minimum cross section S m (t) curve cor- 
responding to the simulation depicted in (a) (the arrow points 
out the time the snapshot was taken), (c) Minimum cross sec- 
tion histogram H(S m ) obtained from tha accumulation of 300 
independent breaking events. 



a given initial size, orientation and temperature) is re- 
flected in its histogram of minimum cross-section H(S m ) 
as a huge peak centred at S m ~ 5 (as it is shown in Fig. 
[TJ;. In those previous works, the relative height of that 
peak or its area has been used as a parameter to classify 
the conditions and probability of formation of icosahedral 
nanowires, as well as the characterization of individual 
pentagonal nanowires was made measuring the length of 
that cross-section plateau. 

However, this measure only provides a qualitative 
value of the average length of the pentagonal nanowires 
formed under stretching, and provides little information 
about the actual lengths of such pentagonal chains. De- 
formations of non-pentagonal regions during the stretch- 
ing or the existence of stages on the formation of the pen- 
tagonal nanowire that provide a cross-section value dis- 
tant from S m = 5 lead to the miscalculation of the length 
of the pentagonal region. This method can not determine 
either the number of pentagonal rings that form the tubu- 
lar structure. In order to overcome these limitations, in 
this paper we present an algorithm that allows the auto- 
matic identification of pentagonal rings structures as well 
as the determination of the actual pentagonal nanotube 
length L p . With this new tool we have revisited the Ni 
case case complementing previous statistical analysis. 



COMPUTATIONAL METHODS 

In our study we have followed the same approach of 
previous papers on Ni to perform the simulated nanowire 
breakage. In short, we have used semi-classical Molecu- 
lar Dynamics (MD) methods at constant temperature to 



study the structure and rupture of metallic nanowires. 
The atomic interactions are described by the Mishin et 
al. parameterization [14] of the embedded atom method 
(EAM) potentials [THIUS]. A nanowire breakage simu- 
lation starts with a parallelepiped of atoms ordered ac- 
cording to a FCC structure with bulk Ni lattice param- 
eter a = 3.52A. The nanowire is stretched along the z 
coordinate at constant velocity till the nanowire breaks. 
During the stretching process the accurate knowledge of 
the atomic coordinates and velocities allows the charac- 
terization of the geometry (in particular the minimum 
cross section S m [T7J [TB]), forces, and kinetic and po- 
tential energies of the stretched nanocontact. The pro- 
cedure of the breakage simulation and the S m calcu- 
lation in atoms units are described in detail in Refs. 

na musing eh is]]. 

We are interested in this study on the appear- 
ance and identification of icosahedral structures, so 
we have focussed on in nanowires stretched along the 
[100] direction due to their large probability to form 
pentagonal nanowires [TUl [HI H2J. The number of 
atoms of the initial nanowire is 1029 (21 layers x 
49 atoms per layer) and the temperatures used were 
T = 4, 160, 225, 300, 375, 465, 550, 690 and 865K (bulk Ni 
melting temperature T m — 1730K). 

In previous works pentagonal staggered nanowires were 
tracked by inspecting the minimum cross-section trace in 
the region S m (t) ~ 5; we define At$ as the time spent by 
the nanowire in the 4 < S m (t) < 6 region (see Fig. [TJd). 
This procedure presents two problems. First, it did not 
allow an exact identification of the simulations creating 
pentagonal structures. Second, as it is discussed in the 
introduction, that parameter, multiply by the stretch- 
ing velocity, could not reflect the actual length of the 
icosahedral nanowire. To overcome these difficulties we 
have developed an algorithm that identifies the pentag- 
onal structures that form the icosahedral nanowire. The 
algorithm also allows us to define its length L p as the 
distance between the outermost pentagonal rings and to 
count the number of pentagonal rings n p that form it. 

The algorithm is based in the determination of the 
angular distribution of the nearest-neighbors atoms and 
provides a parameter (a(z)) which compares the angu- 
lar distribution of the projected nanowire atomic coor- 
dinates with that corresponding to a perfect pentagonal 
nanowire. For a given z coordinate we consider a slice 
perpendicular to the z (stretching) direction with a thick- 
ness of 2 A and centred on the z value. The N t atoms 
inside such slice are projected onto the xy-plane, each one 
getting new 2D coordinates pi; then the centroid of this 
structure is calculated po — ^Ipi/Nt- The angular dis- 
tribution is calculated from the angles ipij between the 
pairs of vectors p/ and p! defining the projected atomic 
coordinates with respect to the centroid (p/ = pi — po). 
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The parameter a is calculated as 

2 ^ \tfij - m<p \ 

where N a is the number of pair of atoms considered, 
ipo = 7r/5 is the reference angle of a perfect staggered 
pentagonal structure, and m is the integer that mini- 
mizes the expression \tpij — mtpo\. To avoid spurious 
contributions from atoms near the centroid, only vectors 
satisfying > a/4 are considered, i.e., center atoms 
are excluded from the calculation of a. This algorithm is 
applied along the z-coordinate of the nanowire, displac- 
ing the imaginary slab Sz — 0.1 A at a time. This results 
in a a(z) pentagonal profile of the nanowire. 

In order to minimize artifacts, the a(z) curve is soft- 
ened over a wider window. This softened curve {ct)(z) is 
defined as 

(a)(z)=— / a(z')dz', (2) 

where a value of Az = lA has been found appropri- 
ate. This average of a ((a)) over a lA interval provides a 
quantity that distinguishes between pentagonal and non- 
pentagonal structures through the nanowire. We have 
observed that if the parameter (a)(z) < 0.5, the set of 
atoms around z forms a structure similar to that of a 
pentagonal ring. On the contrary, if (a)(z) > 0.5 the 
set of atoms presents another structure (bulk like -FCC-, 
helical or or disordered). 

To check the ability of the algorithm to discriminate 
between different structures in Fig. [2] we show the aver- 
age a values obtained using the algorithm along differ- 
ent nanowires with increasing amount of disorder. The 
test structures for the algorithm were: square nanowires 
with atoms sequence 5-4-5-4 taken from a FCC structure 
along the [100] direction; staggered pentagonal nanowires 
with atoms sequence 1-5-1-5; staggered hexagonal 1-6-1- 
6 nanowires; and staggered heptagonal nanowires with 
1-7-1-7 sequence. The degree of disorder was measured 
with the mean displacement a of the atoms from their 
position in the perfect structure: a = J^i Pi ~ P*i,o/RNt) 
where pifi corresponds to the ordered structure atom po- 
sition, and R is the effective radius of the ordered test 
configuration. The initial ordered structures (tr = 0) are 
depicted as inset in Fig. [2] and only in the case of pentag- 
onal nanowire the parameter a takes value 0, being ~1 
for the other structures. As the disorder increases, the a 
average (a) varies: it increases for pentagonal nanowires 
and slightly decreases for the other nanowires. If the dis- 
order with respect the initial structure is strong enough, 
the average of a for all the test nanowires converges to a 
value ~0.9. 

In Fig. [3] we present the result of the algorithm applied 
over a simulated nanowire that presents an icosahedral 
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FIG. 2: The average of the a-parameter (a) versus different 
strengths of the disorder parameter a (see text) for four test 
configurations: square 5-4-5-4, staggered pentagonal 1-5-1-5, 
staggered hexagonal 1-6-1-6 and staggered heptagonal 1-7-1-7 
nanowires (the inset shows the perfect ordered configurations 
of the four test nanowires). a is the mean atomic displace- 
ment of atoms with respect to the perfect position of the 
ordered configuration. The average value of a for disordered 
nanowires was obtained averaging over 50 configurations. 
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FIG. 3: a-parameter (dots) and its average (a) (solid line) 
along the Ni[100] simulated nanowire at T = 375K shown in 
Fig. [l^. For this particular case we include the icosahedral 
nanowire length L p (=21.1 A). The doted line (a) = 5 is the 
limit value we have chosen for the identification of pentagonal 
structures. The minimum values of (a) correspond to the 
position of pentagonal rings forming the icosahedral nanowire. 



structure (that showed in Fig. [Tp. Here we show both 
a(z) and the softened (a)(z) curves. As illustrated in the 
figure, the algorithm returns value near to 1 when is ap- 
plied to the ordered regions of the nanowire, and values 
below 1 for the thinnest part of the nanowire. Minima of 
(a) correspond to the position of the pentagonal rings; as 
they are not perfect ordered structures (but still they can 
be recognized as pentagons) their (a) values are greater 
than cero. We have chosen the value of (a) = 0.5 as the 
limit value to recognized a pentagonal structure. As it 
can be seen in Fig. [2j non-pentagonal tubular structures 
(even with strong disorder) have a (a) value higher than 
0.5. The pentagonal nanowire, even with a relative strong 
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disorder, presents a (a) lower than 0.5; a disordered pen- 
tagon with (a) > 0.5 can not be identified as a regular 
pentagon (see Fig. [2|. The value (a) — 0.5 discriminates 
between pentagonal and non-pentagonal structures. We 
define the pentagonal nanotube length L p (t), observed 
during stretching at a given time t, from the maximum 
and minimum z-coordinates with (a) = 0.5 as it is shown 
in Fig. [3j L™ is the maximum value of L p (t), observed 
when the nanowire is about to break, and n p is the num- 
ber of pentagonal rings forming the icosahedral nanowire 
at its late stage (equivalent to the number of (a) minima 
below 0.5). 

RESULTS 

Once the ability of the algorithm to identify pentago- 
nal structures has been checked, we apply the algorithm 
to our simulations on Ni [100] nanowires. As in previous 
works the aim is to carry out a statistical study for a 
broad range of temperature, obtaining probability distri- 
butions of pentagonal nanowires lengths and number 
of pentagonal rings n p . This will give the optimal tem- 
perature required for maximizing their occurrence prob- 
ability, taking into account that n p = 2 is the minimum 
value defining a pentagonal nanowire. We compare our 
results with those obtained using the time interval Ai 5 
to characterize pentagonal nanowires. 

In Fig. [4] we show the distribution (obtained from 300 
simulations for each temperature) of At 5, L™ and n p for 
a set of temperatures ranging from near to T m /2. 

The At 5 distributions show a smooth shape, being the 
most remarkable differences with the temperature the 
presence of a long tail for the cases ranging from 300 to 
550 K. The initial part of the distribution (up to the max- 
imum position around 30 ps) is due to the plastic defor- 
mation of all the nanowires, and it does not reflect neces- 
sarily the formation of any icosahedral structure. Higher 
values of At 5 would reflect the creation of long structures 
with S m ~ 5, (pentagonal nanowires) and the distribu- 
tion tail would reflect the length of that structures and 
probability of findding one in a stretching proccess. In 
cases with a long tail distribution (300-550 K) the pen- 
tagonal structure in some nanowires seems to last up to 
350 ps, that is equivalent to structures with a length of 
14 A long. 

The way At$ is obtained is similar to that used to 
measure experimental plateau lengths distributions |22j . 
However, we have found that this quantity does not 
reflect the actual maximum length of the icosahedral 
nanowires U p . We will return to this issue later, but 
it is clear when At 5 and L™ distributions are compared. 

The length distributions (Fig. Ep) present a clearly 
peaked shape, with the peaks separated approximately 
by integer values of the distance ds_5 = 2.22A, the cal- 
culated equilibrium separation between successive stag- 



gered pentagonal rings [5] E] ■ The distribution value for 
L p l = corresponds to those cases where the nanowire 
does not form any pentagonal structure (np=0). Those 
cases with < i™ < 2A correspond to nanowires show- 
ing a unique pentagonal ring (n p — 1). We consider 
that the nanowire has formed a pentagonal or icosahe- 
dral nanowire when n p > 2, i.e. when it contains at 
least a full icosahedron in its narrowest part. In agree- 
ment with the At 5 distributions, the L™ distributions 
from T = 300K to 550K show large tails corresponding 
to those pentagonal nanowires including a large number 
of pentagonal rings. In particular, we have detected two 
cases with n p = 11 pentagonal rings for the ensemble of 
breaking events studied at T = 375K set (both pentago- 
nal nanowires were characterized with L™ — 22. 2A and 
23. 7A, respectively). 

In Fig. [IJ; is shown the np distribution for the same 
sets of simulations at different temperatures. Clearly, 
for all the temperatures, the most probable case is that 
the nanowire does not form any pentagonal structure, ie. 
n p = 0. (N(n p ) = agrees with that of 7V(L™ = 0)). 
However for the intermediate range of studied temper- 
atures it is most probable to form a two pentagonal 
rings structure than that only showing a single ring. In- 
deed, for two temperatures (375K and 465K) the most 
probable nanowire structure possesses 3 staggered pen- 
tagons (n p = 3). These two temperatures also exhibit 
the longest icosahedral structures. Ideally every struc- 
ture with a n p > 2 would give a sharp peak in the N^L™) 
distribution. But real structures show tilted pentagonal 
rings, there is not perfect aligned with z-axis and the 
pentagonal rings separation is longer that the equilib- 
rium, producing broad peaks in the N(L™) figure. 

As it has been pointed out previously, the parameter 
Ai 5 does not reflect the actual length of the pentagonal 
nanowire, and its distribution can not been taken as in- 
dication of the actual length distribution. A priori, the 
shortest At$ values must correspond to nonpentagonal 
nanowires. However, for larger At$ values, their max- 
imum length Lp 1 would be shorter than the nanowire 
length increase calculated as Ats times the stretching ve- 
locity. This could be due to the plastic deformation suffer 
by different regions of the nanowire during the stretching 
process. In Fig. [5] we show the L™ value against the cor- 
responding At 5 value for the cases where the stretched 
nanowire creates a pentagonal nanowire (i.e, the plotted 
data are restricted to nanowires with n p > 2). Fig. [5] 
also depicts the expected increased length as function of 
the time (for clarity only results from three temperatures 
are shown; others temperatures shown a similar distribu- 
tion. This figure demonstrates that the former argument 
is wrong; since the maximum icosahedral nanowire length 
is generally larger than that expected froma a direct mea- 
surement of Ate,. 

Fig. [5] proves that, in general, the At$ measure un- 
derestimates the maximum length for long icosahedral 
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FIG. 4: Distribution function for different magnitudes obtained from a set of 300 Ni[100] simulated nanowires with a size 
of 1029 atoms, (a) Time Ats spend in the region of cross section 4.5 < S m (t) < 5.5 by the set of nanowires. (b) Maximum 
pentagonal nanowire length L™; the value N(L™ = 0) corresponds to the number of nanowires from the set of 300 that do not 
show any pentagonal ring in its structure, (c) Number of pentagonal rings n p just before the pentagonal nanowire breakage 
(associates with the length LJT). Each row corresponds to a different temperature: 4, 160, 225, 300, 375, 465, 550, 690 and 865 
K respectively. 



nanowires. For the case of short nanowires Ai 5 tends 
to overestimate the actual pentagonal chain length. The 
magnitude of such underestimation depends on the par- 
ticular nanowire. The cause of the underestimation 
is that, in many nanowires, when the minimum cross- 
section starts to get close to S m ~ 5, there is already 
formed an icosahedal structure with several pentagonal 
rings (generally two or three) (see Fig. [2]d on reference 
|12j for an example). The previous proto-icosahedral 
structure shows S m > 6 and, after stretching, it be- 



comes a classical pentagonal nanowire completely formed 
showing the -1-5-1-5-structure with S m ~ 5. From this 
moment, and upon increasing stretching, the pentago- 
nal nanowire adds new pentagonal rings to the struc- 
ture. Depending on how the formation process occurs 
in every particular event, the At 5 measure lacks to take 
into account the length of several pentagonal rings, and 
therefore underestimated the total icosahedral nanowire 
length. So the pentagonal nanowire length estimations 
of previous studies [10J [TTJ [12] reflected shorter lengths 



6 



5 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
o o 


o 




o 








° o D - ° 

o ° o9 -8°° °°° ° 


^ - 4 m/s 


■ <*» " 

° A n /O ° * 














■ 4 K 

o 375 K 
* 865 K 









50 100 150 200 250 300 350 400 



At 5 (ps) 

FIG. 5: The length versus its corresponding At$ time for 
the Ni[100] simulated nanowires that present the formation 
of an icosahedral nanowire, i.e. n p > 2 (at T = 4, 375 and 
865 K). The solid line is the expected relation between the 
nanowire length increasing and time stretching at a constant 
velocity of 4m/s. 



of icosahedral structures is not a favorable situation for 
low and high temperatures. However, there is a range of 
temperatures (300-550 K) with a large probability (above 
50%) of obtaining long pentagonal nanowires (n p > 2) 
from stretching processes. This temperature dependence 
of the formation of pentagonal chains opens a technolog- 
ical way to optimize the fabrication of these nanoobjects. 
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than the real ones. Only the direct measure of the length 
(applying the algorithm described in this paper or sim- 
ilar methods) can provide information about the forma- 
tion of pentagonal chains, allowing the statistical study 
of the appearance of pentagonal nanowires. 



CONCLUSIONS 

We have present a computational method to identify 
and characterize icosahedral nanowires. This methodol- 
ogy also allows the determination of the pentagonal chain 
length as well as the number of pentagonal rings that 
forms it. We have tested the proposed algorithm for sev- 
eral ad-hoc ordered and disordered structures, proving 
that it can satisfactorily distinguish staggered pentago- 
nal nanowires from other tubular structures. 

The new methodoloy has been applied for statistically 
studying hundreds of Ni[100] breaking nanowire simu- 
lations at different temperatures (ranging from 4K to 
865K), obtaining the pentagonal length and number 
of rings np distributions. We have compared these re- 
sults with the At 5 distributions, already used in previous 
works to identify the existence of pentagonal nanowires. 
We have shown that the quantity At$ generally under- 
estimated the length of the icosahedral nanowire and, as 
consequence it is not adequate for the characterization of 
icosahedral nanowires. 

The Li™ distributions show a peaked structure, similar 
to those found for linear atomic chains [22 with clear 
peaks separated by integer values of the calculated equi- 
librium distance between two consecutive staggered pen- 
tagons. The 



rip distribution shows that the formation 
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